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Preface 


The  Long  Term  Resource  Monitoring  Program  (LTRMP)  was  authorized  under  the  Water  Resources 
Development  Act  of  1986  (Public  Law  99-662)  as  an  element  of  the  Environmental  Management 
Program  for  the  Upper  Mississippi  River  System.  The  LTRMP  is  implemented  by  the  Upper  Midwest 
Environmental  Sciences  Center  of  the  U.S.  Geological  Survey,  in  cooperation  with  the  five  Upper 
Mississippi  River  System  states  (Illinois,  Iowa,  Minnesota,  Missouri,  and  Wisconsin),  with  guidance 
and  Program  responsibility  provided  by  the  U.S.  Army  Corps  of  Engineers. 

The  mission  of  the  LTRMP  is  to  provide  decision  makers  with  information  to  maintain  the  Upper 
Mississippi  River  System  as  a  viable  large  river  ecosystem  given  its  multiple-use  character.  The 
longterm  goals  of  the  Program  are  to  understand  the  system,  determine  resource  trends  and  impacts, 
develop  management  alternatives,  manage  information,  and  develop  useful  products. 
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Abstract:  This  study  was  designed  to  assess  whether  fish  community  data  collected  for  the  Long  Term 
Resource  Monitoring  Program  (LTRMP)  from  six  regional  trend  analysis  (RTA)  areas  of  the  Upper 
Mississippi  River  System  (UMRS)  can  be  used  to  make  inferences  to  the  system  as  a  whole.  Spatial 
coverage  of  fish  monitoring  for  three  LTRMP  field  stations  was  extended  to  "outpools"  immediately 
above  and  below  RTA  pools  4  and  13  and  the  Open  River  Reach  from  June  15  to  October  31, 2000. 
Also,  we  sampled  Navigation  Pools  19  and  20  using  LTRMP  electrofishing  methodology  in  September 
2000.  Multivariate  statistical  analyses  were  used  to  group  pools  on  the  basis  of  fish  community 
composition  and  community  structure.  Cluster  analysis  of  community  composition  and  structure 
data  revealed  two  major  groups  of  pools;  upper  pools  (i.e.,  northern)  and  lower  pools  (i.e.,  southern). 
Navigation  Pools  19  and  20  grouped  with  lower  pools  in  terms  of  community  composition,  and  with 
upper  pools  in  terms  of  community  structure.  Analysis  of  community  composition  data  yielded  four 
subgroups,  with  La  Grange  Pool  forming  its  own  subgroup.  Analysis  of  community  structure  yielded 
five  subgroups,  with  La  Grange  Pool  and  Pool  8  forming  unique  subgroups.  In  general,  all  outpools 
grouped  with  the  nearest  RTA  pools  for  both  community  composition  (no  exception)  and  community 
structure  (one  exception).  Strong  correlations  between  the  community  composition  and  structure 
matrices  with  distance  between  pools  suggest  that  fish  communities  in  relatively  close  pools  are  more 
similar  than  in  pools  separated  by  larger  distances.  Habitat  variables  measured  during  electrofishing 
collections  were  significantly  correlated  with  spatial  variation  of  fish  composition  and  community 
structure,  but  provided  only  marginal  improvements  to  correlations  with  distance  between  pools  alone. 


Results  of  this  study  lend  support  to  the  premise  that  LTRMP  fish  community  data  could  potentially  be 
used  to  make  inferences  to  the  entire  UMRS,  because  current  RTA  areas  are  evenly  distributed  within 
the  major  pool  groupings  identified  in  this  study.  Nevertheless,  further  research  is  needed  to  resolve 
how  fish  communities  in  Navigation  Pools  19  and  20  and  other  lower  UMRS  pools  compare  to  present 
RTA  areas. 

Key  words:  Analysis,  fish  community,  LTRMP,  Mississippi  River,  navigation  pool,  trend  analysis 
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Introduction 

The  Long  Term  Resource  Monitoring  Program 
(LTRMP)  was  authorized  by  the  Water  Resources 
Development  Acts  of  1986  and  1999  as  an 
element  of  the  U.S.  Army  Corps  of  Engineers' 
Environmental  Management  Program.  The  primary 
mission  of  the  LTRMP  is  to  provide  resource 
managers  with  the  information  needed  to  maintain 
the  Upper  Mississippi  River  System  (UMRS)  as 
a  viable  multiple-use  ecosystem.  Four  long-term 
goals  established  for  the  LTRMP  are  (1)  increasing 
our  understanding  of  how  the  river  ecosystem 
operates.  (2)  monitoring  UMRS  natural  resources 
status  and  trends,  (3)  assisting  in  the  evaluation  of 
management  alternatives,  and  (4)  managing  and 
providing  access  to  resulting  data,  information,  and 
products  (U.S.  Army  Corps  of  Engineers  1997). 
Standardized  monitoring  of  water  quality,  aquatic 
vegetation,  aquatic  invertebrates,  and  fishes  are  key 
components  of  the  LTRMP. 

The  vast  geographic  expanse  of  the  UMRS, 
which  includes  1,300  miles  of  navigable  rivers 
and  a  basin  encompassing  190,000  square  miles 
including  numerous  tributaries,  presents  a  major 
challenge  to  the  LTRMP.  The  design  of  any 
monitoring  program  focused  on  a  large  ecosystem 
must  find  a  balance  between  covering  enough 
spatial  area  to  allow  for  system  wide  inferences  and 
the  ability  to  obtain  enough  detailed  information 
to  describe  and  understand  the  interaction  of 
ecosystem  components.  This  challenge  is  even 
greater  when  populations  and  communities  of 
organisms  are  a  primary  focus  of  a  monitoring 
program,  given  that  population  dynamics  of 
different  species  within  an  ecosystem  usually 
operate  at  several  different  spatial  and  temporal 
scales  (Wiens  1989;  Levin  1992).  Therefore, 
monitoring  programs  designed  to  track  population 


variation  of  several  species  need  to  sample  multiple 
spatial  subunits  appropriate  for  populations 
operating  at  small  scales,  with  sufficient  replication 
of  subunits  to  track  populations  operating  at  larger 
scales  and  allow  inferences  to  the  entire  system. 

For  a  system  as  spatially  extensive  as  the 
UMRS,  it  is  clear  that  many  species  have  multiple 
populations  within  this  river-floodplain  system. 
Because  a  major  goal  of  LTRMP  is  to  provide 
system  wide  inferences  for  the  monitored 
components,  it  is  also  clear  that  the  monitoring 
design  must  include  substantial  spatial  coverage. 
The  original  plans  for  the  LTRMP  included  sample 
collection  from  22  river  reaches  (U.S.  Army 
Corps  of  Engineers  1997).  Logistical  constraints, 
however,  reduced  the  number  of  areas  sampled 
from  22  river  reaches  to  6  regional  trend  analysis 
areas  (RTA;  Navigation  Pools  4,  8,  13,  and  26 
and  the  Open  River  Reach  on  the  Mississippi 
River  and  La  Grange  Pool  on  the  Illinois  River). 
Lubinski  et  al.  (2001)  conducted  a  power  analysis 
demonstrating  that  the  present  LTRMP  design 
provides  adequate  statistical  power  to  detect 
inter-annual  variation  for  most  water  quality 
parameters  and  many  fishes  within  the  six  RTA 
areas.  Nevertheless,  how  well  these  six  RTA  areas 
reflect  the  overall  biotic  and  abiotic  conditions  of 
the  entire  UMRS  remains  unknown. 

This  study  was  initiated  to  expand  the 
spatial  coverage  of  LTRMP  fisheries  monitoring 
(specifically,  electrofishing  and  seining)  to  outpools 
immediately  above  and  below  three  RTA  areas  in 
an  attempt  to  provide  further  insight  into  the  ability 
of  LTRMP  to  make  ecological  inferences  to  the 
UMRS  as  a  whole.  We  present  statistical  analysis 
of  this  study,  examining  spatial  variation  of 
community  composition  (the  presence  or  absence 
of  species)  and  community  structure  (relative 
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abundance  of  species)  of  UMRS  fishes.  Four  major 
goals  are  addressed: 

1 .  Determine  how  sampled  pools  group  based 
on  community  composition  of  fishes. 

2.  Determine  how  sampled  pools  group  based 
on  community  structure  of  fishes. 

3 .  Determine  where  outpools  group  relative  to 
adjacent  RTA  areas. 

4.  Determine  whether  groupings  of  pools 
based  on  community  composition  and 
structure  correspond  to  spatial  variation  of 
important  habitat  factors. 

Methods 

Fish  Sampling 

This  study  extended  the  spatial  coverage  of 
fish  monitoring  for  three  LTRMP  monitoring 


Gutreuter  et  al.  (1995)  described  standard 
LTRMP  methodology  for  electrofishing  and 
seining  in  detail.  Sampling  locations  were 
selected  using  a  stratified  (by  habitat  type)  random 
design  (Lubinski  et  al.  2001).  Electrofishing  was 
conducted  using  pulsed-DC  output  with  two  ring 
anodes  and  the  boat  hull  serving  as  the  cathode 
and  voltage  and  amperage  were  adjusted  for  water 
temperature  and  conductivity  to  achieve  a  power 
output  of  3,000  W.  Two  dippers  collected  fish. 
Electrofishing  was  conducted  along  shorelines 
continuously  for  15  min  at  each  sample  collection 
site.  Data  on  water  temperature,  depth  (average 
for  each  collection  site),  conductivity,  and  habitat 
were  collected  with  each  sample  (Table  1).  Seining 
was  conducted  using  a  10.7-m-long  3-mm-mesh 
bag  seine.  Seines  are  fished  along  banks  in  water 
<1.2  m.  One  end  of  the  seine  was  anchored  to  the 
bank  and  the  other  end  was  deployed  perpendicular 
to  the  bank  and  swept  downstream.  All  fish  were 


locations.  Mississippi  River  navigation 
pools  inunediately  above  and  below 
RTA  pools  4  and  13  and  Open  River 
Reach  (Figure  1)  were  sampled  using 
standard  LTRMP  electrofishing 
methodology  from  June  15  to  October 
31, 2000.  Throughout  this  report,  these 
areas  will  be  referred  to  collectively  as 
outpools  and  individually  as  Pools  3, 

5, 12, 14, 29,  and  31  (note;  29  and  31 
are  nonpooled  river  reaches).  Standard 
LTRMP  monitoring  also  occurred 
in  all  six  RTA  pools  (4,  8, 13, 26, 

La  Grange,  and  Open  River  Reach) 
during  2000  and  electrofishing  samples 
were  collected  from  Pools  19  and  20 
by  all  LTRMP  fisheries  personnel  on 
September  12  and  13, 2000  (Figure  1). 
We  were  interested  in  examining  how 
fish  communities  in  Pools  19  and 
20  grouped  with  other  UMRS  pools 
because  Lock  and  Dam  19  is  believed  to 
present  a  barrier  to  migrations  of  certain 
fishes  (Kelner  and  Seitman  2000).  Plans 
for  this  study  called  for  seining  in  all 
outpools.  Unfortunately,  this  proved  to 
be  logistically  impractical  and  only  was 
accomplished  for  two  outpools  (12  and 

14). 


Figure  1.  Map  of  the  Upper  Mississippi  River  System  showing  the  six  regional  trend 
analysis  pools  monitored  by  the  Long  Term  Resource  Monitoring  Program  and  the  eight 
outpools  sampled  with  electrofishing  during  2000. 
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Table  1.  Habitat  variables  routinely  collected  from  each  electrofishing  site  for  the  Long  Term  Resource  Monitoring 
Program  (Gutreuter  et  al.  1995), 


Habitat  factor 

Units 

Explanation 

Seech i 

cm 

Measurement  of  water  transparency  in  cm 

Conductivity 

S/cm 

Conductivity  measured  to  the  nearest  I  S/cm 

Flow 

m/scc 

Rate  at  which  the  water  is  flowing  given  in  m/sec 

Temperature 

C 

Temperature  of  the  water  in  C 

Depth 

m 

Water  depth  in  fractions  of  meters 

Emergent/Submersed  vegetation 

0,  1,2,3 

0  =  0%  coverage;  1  =  1-19%  coverage; 

2  =  20-49%  coverage;  3  =  50%  coverage 

Vegetation  density 

0, 1,2 

0  =  no  veg;  1  -  sparse;  2  =  dense 

Substrate 

1,2,  3,4 

1  =  silt;  2  =  silt/clay/little  sand;  3  =  sand/mostly  sand; 

4  =  gravel/rock/hard  clay 

Woody  structure 

pres/abs 

presence  or  absence  of  woody  structure 

Revetment 

pres/abs 

presence  or  absence  of  shoreline  revetment 

Inlet/Outlet 

pres/abs 

presence  or  absence  on  an  inlet/outlet  channel  to  a  backwater  lake 

Flooded  terrestrial  vegetation 

pres/abs 

presence  or  absence  of  flooded  terrestrial  vegetation 

identified,  measured,  and  enumerated  following 
standard  LTRMP  protocol  (Gutreuter  et  al.  1995). 

Statistical  Analysis 

We  examined  spatial  variation  in  fish 
community  composition  and  structure  among  the 
six  RTA  pools  and  eight  outpools.  Community 
composition  refers  to  the  presence  or  absence  of 
species,  whereas  community  structure  refers  to  the 
abundance  of  species  as  measured  by  mean  catch- 
per-unit-effort  (CPUE  equals  number  per  15  min, 
weighted  by  habitat  strata).  Separate  analyses  were 
conducted  for  electrofishing  and  seining  data, 
and  all  analyses  were  conducted  using  SAS  for 
Windows  (SAS  Institute,  Inc.  1999)  and  Primer  for 
Windows  (Primer-E  LTD  2001).  Our  analysis  of 
seining  data  was  limited  to  community  composition 
because  the  power  to  detect  variation  in  abundance 
of  fishes  from  LTRMP  seining  data  differs  greatly 
among  pools  (Lubinski  et  al.  2001). 

For  both  response  variables  (presence/absence, 
CPUE),  we  used  cluster  analysis  and  nonmetric 
multidimensional  scaling  (NMDS)  to  identify 
groupings  of  pools.  These  analyses  were  based 
on  a  Euclidian  distance  matrix  for  community 
composition  data,  and  a  Bray-Curtis  similarity 


matrix  for  community  structure  data.  Catch-per- 
unit-effort  data  were  square-root  transformed 
to  better  conform  to  multivariate  normality 
assumptions.  This  transformation  also  dampens  the 
influence  of  very  abundant  species  for  conrununity 
structure  analysis  (Clarke  and  Warwick  1994). 

We  limited  the  community  structure  analysis  to 
16  species  for  which  electrofishing  had  power 
0.80  to  detect  a  20%  interannual  abundance 
change  in  at  least  one  habitat  strata  of  an  RTA  pool 
based  on  the  Lubinski  et  al.  (2001)  power  analysis 
of  LTRMP  components  (Table  2).  This  somewhat 
conservative  criterion  was  adopted  to  help  ensure 
that  the  patterns  of  relative  abundance  used  in  these 
analyses  reflect  true  ecological  patterns  rather  than 
sampling  artifacts.  Hybrids  and  fish  not  identified 
to  species  were  omitted  from  all  analyses. 

Three  criteria  were  used  to  determine  the 
subgrouping  level  in  our  cluster  analysis.  First,  we 
used  rarefaction  curves  from  the  six  RTA  pools 
to  visually  determine  the  minimum  number  of 
individuals  needed  to  reach  the  asymptote  of  the 
rarefaction  curve  (i.e.,  the  sampling  effort  needed 
to  adequately  describe  species  composition). 
Acceptable  subgrouping  levels  should  not  isolate 
undersampled  pools  because  this  isolation 
could  have  resulted  from  a  sampling  artifact. 
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Table  2.  Mean  abundance  (square  root  #/ 1 5  min)  of  the  1 6  species  used  in  community  structure  analysis^ 


Species" 

Upper  Mississippi  River  System  Navigation  Pool 

3 

4 

5 

8 

12 

13 

14 

19 

20 

LG 

26 

29 

OR 

31 

Gizzard  shad 

7.50 

6.30 

3.53 

2.18 

3.45 

4.12 

2.94 

2.61 

2.74 

6.98 

4.51 

3.23 

5.36 

4.14 

Emerald  shiner 

7.13 

3.99 

1.77 

2.34 

3.93 

3.99 

1.69 

6.85 

6.96 

1.08 

1.49 

1.72 

2.46 

1.26 

Common  carp 

2.07 

2.23 

2.19 

1.18 

2.23 

2.30 

2.30 

2.11 

1.79 

2.76 

2.46 

1.99 

1.20 

1.71 

Bluegill 

0.52 

2.97 

3.09 

6.28 

3.19 

3.00 

3.76 

1.48 

0.17 

1.75 

0.84 

0.05 

0.09 

0.08 

Freshwater  drum 

1.44 

0.67 

0.56 

0.34 

0.82 

0.97 

0.91 

1.86 

1.46 

1.50 

0.95 

1.26 

1.22 

0.97 

Largemouth  bass 

0.05 

1.43 

1.21 

2.66 

1.88 

2.27 

2.24 

1.22 

0.12 

1.23 

0.26 

0.00 

0.04 

0.02 

Spotfin  shiner 

1.91 

1.06 

1.12 

2.62 

1,18 

1.14 

0.21 

1.21 

1.53 

0.00 

0.10 

0.08 

0.00 

0.00 

Bullhead  minnow 

1.16 

0.75 

0.89 

3.08 

2.03 

1.46 

0.89 

0.81 

0.36 

0.11 

0.13 

0.00 

0.01 

0.02 

White  bass 

1.45 

0.63 

0.39 

0.45 

0.76 

0.76 

0.41 

0.89 

0.72 

1.63 

0.49 

0.69 

0.72 

0.46 

Channel  catfish 

0.26 

0.09 

0.12 

0.09 

0.58 

0,69 

0.38 

1.77 

0.55 

1.17 

0.71 

0.74 

0.83 

0.69 

Black  crappie 

0.08 

0.92 

0.90 

0.95 

0.43 

1.14 

1.01 

0.10 

0.08 

0.37 

0.06 

0.02 

0.00 

0.00 

Smallmouth  buffalo 

0.22 

0.20 

0.15 

0.23 

0,84 

0.24 

0.12 

0.22 

0.63 

1.68 

0.59 

0.08 

0.42 

0.40 

Shorthead  redhorse 

0.82 

1.16 

1.18 

0.91 

0.55 

0.28 

0.12 

0.13 

0.07 

0.12 

0.06 

0.08 

0.04 

0.00 

Smallmouth  bass 

0.92 

0.73 

0.98 

0.65 

0.46 

0.06 

0.18 

0.10 

0.05 

0.00 

0.00 

0.01 

0.02 

0.00 

Silver  redhorse 

0.14 

0.84 

1.09 

0.84 

0.00 

0.00 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Bigmouth  buffalo 

0.02 

0.10 

0,19 

0.00 

0.42 

0.41 

0.15 

0.13 

0.00 

1.05 

0.16 

0.01 

0.17 

0.10 

•Samples  were  collected  from  June  15  through  October  31, 2000.  LG  =  La  Grange  Pool  of  the  Illinois  River.  OR  =  Open  River 
Reach.  Note:  Pools  29  and  31  are  open  river  reaches. 

•Gommon  names  for  fishes  follow  Robins  et  al.  (1991). 


Second,  we  calculated  the  mean  and  range  of 
Euclidian  distance  (community  composition) 
and  Bray-Curtis  similarity  (community  stmcture) 
between  consecutive  years  for  each  RTA  pool 
using  LTRMP  data  from  1994  to  2000.  Many 
of  the  differences  between  consecutive  years  in 
community  composition  and  structure  data  from 
LTRMP  electrofishing  can  be  attributed  to  sampling 
artifacts  (i.e.,  electrofishing  does  not  sample  all  the 
species  present  in  a  pool  each  year).  Therefore,  the 
critical  Euclidean  distance  for  deriving  robust  pool 
groupings  should  be  greater  than  the  range  of  these 
year-to-year  differences.  Finally,  we  accepted  only 
subgrouping  levels  that  produced  groups  that  were 
easily  illustrated  using  NMDS  in  either  two  or  three 
dimensions  with  a  stress  value  <  0.05.  Stress  value 
is  a  measure  of  "goodness-of-fit"  for  NMDS  with 
small  values  indicating  a  better  fit  than  large  stress 
values  (Clarke  and  Warwick  1994).  Because  seining 
data  were  available  only  for  a  small  subset  of  pools, 
we  determined  subgrouping  using  NMDS  criteria 
alone. 


Because  cluster  analysis  and  NMDS  are 
data  exploration  techniques,  we  used  analysis 
of  similarity  (ANOSM)  to  test  for  significant 
variation  in  fish  community  composition  and 
structure  among  groups.  Analysis  of  similarity 
is  analogous  to  univariate  ANOVA  in  that  it  tests 
for  significant  differences  among  groups.  Unlike 
ANOVA,  however,  ANOSIM  uses  Monte  Carlo 
simulation  to  determine  significance  rather  than 
probability  inferences  from  an  assumed  statistical 
distribution.  Nonetheless,  our  use  of  ANOSIM 
here  is  clearly  a  post-hoc  test  and  results  should  be 
interpreted  with  caution.  Also,  we  used  similarity 
breakdown  analysis  (Clarke  and  Warwick  1994; 
SIMPER  procedure  in  Primer-E  LTD  2001)  to 
determine  the  contribution  of  species  to  Bray- 
Curtis  dissimilarity  between  community  structure 
groupings. 

Finally,  we  used  the  electrofishing  data  to 
determine  whether  variation  in  fish  community 
composition  and  structure  among  pools 
corresponded  with  variation  in  habitat  factors.  Two 
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sources  of  habitat  data  available  for  comparisons 
were  those  measured  directly  in  the  field  with  each 
electrofishing  collection  (Table  1)  and  a  suite  of 
aquatic  or  geomorphic  variables  quantified  for 
the  habitat  needs  assessment  (HNA)  query  tool 
that  were  available  for  a  subset  of  pools  (Table  3; 
DeHaan  et  al.  2000;  Koel  2001).  For  both  sets  of 
habitat  variables,  normalized  (mean  =  0,  standard 
deviation  =  1 )  Euclidean  distance  matrices  were 
calculated  and  Mantel  tests  were  used  to  determine 
correlations  with  both  the  Euclidian  matrix  from 
fish  community  composition  data  and  the  Bray- 
Curtis  similarity  matrix  from  community  structure 
data.  A  canonical  Mantel  test  (Clarke  and  Warwick 
1994;  BioEnv  procedure  in  Primer-E  LTD  2001) 
was  used  to  determine  the  combination  of  habitat 
variables  that  would  provide  the  greatest  correlation 
with  community  data.  Because  many  habitat 
variables  can  co-vary  with  latitude,  we  included 
distance  in  river  miles  (relative  to  Navigation 
Pool  3)  in  both  habitat  data  sets  to  help  determine 
whether  correlations  with  habitat  reflect  covariation 
with  the  distance  between  pools.  For  La  Grange 
Pool,  distance  was  the  sum  of  river  miles  between 
Pools  3  and  26,  and  river  miles  between  Pool  26 
and  La  Grange  Pool.  We  also  used  a  Mantel  test 


to  examine  whether  habitat  similarity  (normalized 
Euclidian  distance  matrix  for  all  habitat  features 
measured  in  the  field)  was  correlated  with  distance 
between  pools. 

Results 

Electrofishing 

A  total  of  1 18,139  fishes  were  collected 
comprising  100  species  (Table  4).  The  species 
with  the  greatest  overall  abundance  were  gizzard 
shad,  emerald  shiner,  bluegill,  and  common  carp, 
which  together  accounted  for  over  71  %  of  the  total 
catch.  The  mean  number  of  sj)ecies  captured  in  a 
pool  was  49,  ranging  from  38  to  60.  Rarefaction 
curves  suggest  that  sampling  sufficient  to  collect 
at  least  5,000  fishes  is  needed  to  do  an  adequate 
job  of  describing  community  composition  within 
a  given  pool  or  reach  (Figure  2).  Because  fewer 
than  5,(X)0  fishes  were  collected  from  several  of  the 
outpools  (Figure  3),  we  rejected  any  subgrouping 
level  that  isolated  these  undersampled  pools  in  our 
cluster  analysis  of  community  composition  and 
community  structure  data. 


Table  3.  Surface  area  (hectares)  of  aquatic/geomorphic  habitat  variables  as  defined  by  the  habitat  needs  assessment  query 
tool  for  Upper  Mississippi  River  Navigation  Pools  4-26  and  open  river  reaches  (29,  OR,  31).* 


Area  (ha)  of  aquatic/geomorphic  habitat  type 

Pool 

MNC 

MCB 

TWZ 

SCH 

TCH 

TRC 

CFL 

CFS 

CIM 

TIS 

CTF 

TOC 

4 

1,093 

448 

12 

463 

2 

97 

10,320 

1,567 

408 

1,848 

8,438 

24,695 

5 

337 

536 

22 

278 

0 

59 

155 

733 

2,178 

863 

6,025 

11,187 

8 

627 

603 

21 

510 

1 

30 

1,125 

1,573 

4,024 

2,966 

3,478 

14,957 

12 

596 

1,506 

26 

740 

20 

4 

401 

545 

864 

1,446 

1,645 

7,794 

13 

1,569 

1,141 

20 

789 

105 

32 

1,242 

1,902 

3,556 

2,414 

8,494 

21,262 

14 

561 

2,127 

22 

599 

4 

24 

668 

0 

0 

1,357 

3,107 

8,470 

19 

1,350 

5,273 

30 

1,527 

1 

93 

868 

1,282 

1,069 

2,297 

14,033 

27,823 

20 

574 

1,728 

39 

545 

4 

200 

23 

0 

0 

786 

4,829 

8,727 

26 

1,467 

2,875 

28 

1,483 

14 

51 

409 

0 

245 

2,530 

18,663 

27,764 

•Habitat  variables  were  the  main  navigation  channel  (MNC),  main  channel  border  (MCB),  tailwater  (TWZ), 
secondary  channel  (SCH).  tertiary  channel  (TCH),  tributary  channel  (TRC),  contiguous  floodplain  lake  (CFL), 
contiguous  floodplain  shallow  aquatic  area  (CFS),  contiguous  impounded  area  (CIM),  terrestrial  island  (TIS), 
contiguous  terrestrial  floodplain  (CTF),  and  total  contiguous  habitat  area  (TOC).  Data  were  not  available  for 
Pool  3,  La  Grange  Pool  of  the  Illinois  River,  or  open  river  reaches  29  and  31 . 
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Upper  Mississippi  River  System  Navigation  Pools^ 
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Pumpkinseed 


Upper  Mississippi  River  System  Navigation  Pools'* 
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Slenderhead  darter 


Upper  Mississippi  River  System  Navigation  Pools^ 
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Blackspotted  topminnow 


Upper  Mississippi  River  System  Navigation  Pools^ 
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Number  of  Individuals  Sampled 

Figure  Z  Rarefaction  curves  from  electrofishing  data  for  the  six  regional  trend  analysis  pools  in  2000.  At  least  5,000 
individuals,  it  seems,  should  be  sampled  to  adequately  describe  fish  community  composition  (i.e.,  to  reach  the  rarefaction 
curve  asymptote). 


Figure  3.  The  total  number  of  fish  captured  by  electrofishing  from  the  six  regional  trend  analysis  pools  and  eight  outpools 
during  2000.  The  horizontal  line  depicts  the  minimum  number  of  individuals  (5,000)  needed  to  be  sampled  to  adequately 
describe  community  composition  (see  Figure  2). 
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Cluster  analysis  of  community  composition  As  with  community  composition,  cluster 

data  revealed  two  major  groupings  of  pools:  upper  analysis  based  on  community  structure  of  fishes 

and  lower  pools  (Figure  4A).  Twenty-eight  species  revealed  two  major  subgroups  of  pools,  upper 

were  captured  only  within  lower  pools  and  18  only  and  lower.  In  this  analysis,  however,  Pools  19  and 

within  upper  pools.  The  average  Euclidian  distance  20  grouped  with  the  upper  pools  rather  than  the 

between  consecutive  years  at  an  RTA  pool  was  lower  pools  (Figure  5A).  Six  species  accounted 
3.60,  and  ranged  from  3.16  to  4.47.  Therefore,  we  for  more  than  70%  of  the  dissimilarity  between 

chose  5.0  as  our  subgrouping  distance,  resulting  upper  and  lower  pool  groupings.  Upper  pools  were 

in  four  subgroups:  all  upper  pools.  La  Grange  characterized  by  greater  abundance  of  emerald 

Pool,  the  Open  River  Reach  and  Pools  29  and  shiner,  bluegill,  largemouth  bass,  bullhead  minnow, 

31,  and  Pools  19,  20,  and  26  (Figure  4A).  These  and  spotfin  shiner  relative  to  lower  pools.  Lower 

four  groups  were  illustrated  by  NMDS  in  three  pools  had  greater  abundance  of  gizzard  shad 

dimensions  with  a  stress  value  =  0.02  (Figure  4B).  compared  with  upper  pools. 

Analysis  of  similarity  revealed  significant 

differences  between  the  two  major  groupings  The  mean  Bray-Curtis  similarity  value  for 

(Global  R  =  0.83;  P  =  0.001),  and  among  the  four  consecutive  years  in  RTA  pools  was  85.7  (range 

subgroups  (Global  R  =  0.95;  P  =  0.001).  from  74.4  to  92.6;  Figure  5A).  We  chose  a 


Figure  4.  Cluster  analysis 
of  fish  community  composition 
data  collected  through 
electrofishing.  The  lower  solid 
line  and  two  dashed  lines 
indicate  the  mean  and  range 
of  Euclidean  distance  for 
consecutive  years  in  regional 
trend  analysis  pools.  The 
upper  solid  line  indicates  the 
subgrouping  level  (Euclidean 
distance  of  5)  resulting  in  four 
subgroups  (circled). 

/'Sj  Nonmetric  multidimensional 
scaling  plot  of  Upper 
Mississippi  River  System  in 
three  dimensions.  The  four 
subgroups  from  the  cluster 
analysis  are  circled. 
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Figure  5.  Cluster  analysis 
of  fish  community  structure 
data  collected  through 
electrofishing.  The  lower  solid 
line  and  two  dashed  lines 
indicate  the  mean  and  range 
of  Bray-Curtis  similarity  for 
consecutive  years  in  regional 
trend  analysis  pools.  The 
upper  solid  line  indicates  the 
subgrouping  level  (Bray-Curtis 
similarity  =  70)  resulting  in  five 
subgroups  (circled). 

(Bj  Nonmetric  multidimensional 
scaling  plot  of  Upper 
Mississippi  River  System  in 
three  dimensions.  The  five 
subgroups  from  the  cluster 
analysis  are  circled. 


The  mean  Bray-Curtis  similarity  value  for 
consecutive  years  in  RTA  pools  was  85.7  (range 
from  74.4  to  92.6;  Figure  5/4).  We  chose  a 
Bray-Curtis  value  of  70  for  our  subgrouping 
distance,  producing  five  subgroups:  Pools  26, 
29,  and  31  and  the  Open  River  Reach  (subgroup 
A);  La  Grange  Pool  (subgroup  B);  Pools  3,  19, 
and  20  (subgroup  C);  Pools  4,  5, 12, 13,  and 
14  (subgroup  D);  and  Pool  8  (subgroup  E). 
Nonmetric  multidimensional  scaling  illustrates 
these  five  subgroups  in  three  dimensions  with  a 
stress  value  =  0.03  (Figure  5B).  Eleven  species 
contributed  more  than  70%  to  the  dissimilarity 
among  subgroups  (Table  5).  Pool  8  (subgroup 
E)  had  the  greatest  abundance  of  black  crappie. 


bullhead  minnow,  bluegill,  largemouth  bass,  and 
spotfin  shiner.  La  Grange  Pool  (subgroup  B) 
had  the  greatest  abundance  of  bigmouth  buffalo, 
common  carp,  smallmouth  buffalo,  and  white  bass. 
Subgroup  C  had  the  greatest  abundance  of  emerald 
shiner  and  freshwater  drum.  Analysis  of  similarity 
revealed  significant  differences  between  upper  and 
lower  pool  groupings  (Global  R  =  0.67;  P  =  0.001 ) 
and  among  the  five  subgroups  (Global  R  =  99; 

P  =  0.001). 

Habitat  Correlations 

Strong  correlations  between  community 
composition  (r  =  0.80)  and  structure  (r  =  0.63) 
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Table  5.  Mean  abundance  (square  root  #/ 1 5  min)  of  the  1 1  species  that  contributed  more  than  70%  to  the  dissimilarity  among  pool 
subgroups  identified  through  cluster  and  non-metric  multidimensional  scaling  analysis  of  eiectrofishing  community  structure  data. 


Subgroups  and  UMRS  Pools 

A 

B 

C 

D 

E 

Species" 

26-31 

La  Grange 

3,19,20 

4,5,12,13,14 

8 

Black  crappie 

0.02 

•  0.37 

0.09 

0.88 

0.95- 

Bluegill 

0.26 

1.75 

0.72 

3.20 

6.28 

Bigmouth  buffalo 

0.11 

1.05 

0.05 

0.25 

0.00 

Bullhead  minnow 

0.04 

0.11 

0.78 

1.20 

3.08 

Common  carp 

1.84 

2.76 

1.99 

2.25 

1.18 

Emerald  shiner 

1.73 

1.08 

6.98 

3.07 

2.34 

Freshwater  drum 

1.10 

1.50 

1.59 

0.79 

0.34 

Gizzard  shad 

4.31 

6.98 

4.28 

4.07 

2.18 

Smallmouth  buffalo 

0.37 

1.68 

0,35 

0.31 

0.23 

Spotfin  shiner 

0.04 

0.00 

0.94 

0.94 

2.62 

White  bass 

0.59 

1.63 

1.02 

0.59 

0.45 

•The  maximum  abundance  for  each  species  are  in  bold. 
'’Common  names  for  fishes  follow  Robins  et  al.  (1991). 


matrices  with  distance  between  pools  suggest 
that  fish  communities  in  relatively  close  pools 
are  more  similar  than  pools  separated  by  larger 
distances  (Figures  6A-6B).  Also,  habitat  similarity 
was  correlated  (r  =  0.74)  with  distance  between 
pools  (Figure  6Q  and  most  of  the  habitat  variables 
measured  during  electrofishing  sampling  showed 
longitudinal  variation  (Figure  7).  These  habitat 
variables  were  significantly  correlated  with  both 
community  composition  (r  =  0.75;  P  =  0.0001) 
and  community  structure  (r  =  0.64;  P  =  0.0010). 
The  canonical  Mantel  procedure  revealed 
that  the  strongest  correlations  (r  =  0.90)  with 
community  composition  were  with  a  habitat  matrix 
composed  of  distance  between  pools,  water  depth, 
conductivity,  vegetation  density,  and  the  frequency 
of  woody  stmcture  presence.  The  strongest 
correlations  (r  =  0.73)  for  community  stmcture 
were  with  a  habitat  matrix  composed  of  distance 
between  pools,  flow,  vegetation  density,  frequency 
of  woody  stmcture  presence,  and  the  frequency  of 
flooded  terrestrial  vegetation  presence.  Note  that 
these  correlations  are  only  marginal  improvements 
over  correlation  with  distance  between  pools  alone. 

Mantel  tests  revealed  significant  correlations 
between  HNA  variables  with  both  community 
composition  (r  =  0.46;  P  =  0.0083)  and  community 
stmcture  (r  =  0.37;  P  =  0.0225).  Nevertheless, 


correlations  with  distance  between  pools  alone 
were  stronger  (r  =  0.79;  r  =  0.69)  and  the  canonical 
Mantel  procedure  was  unable  to  add  any  HNA 
habitat  variables  that  could  more  than  trivially 
improve  these  correlations. 

Seining 

A  total  1 15,820  fishes  from  81  species  were 
captured.  The  species  with  the  greatest  overall 
abundance  were  emerald  shiner,  mimic  shiner, 
river  shiner,  bluegill,  bullhead  minnow,  and  gizzard 
shad.  Together,  these  six  species  accounted  for 
over  80%  of  the  total  catch.  Cluster  analysis  of 
seining  data  revealed  two  major  groupings  of 
pools:  upper  (4,  8, 12, 13,  and  14)  and  lower  pools 
(26,  La  Grange,  and  Open  River  Reach;  Figure  8A). 
These  two  groupings  were  illustrated  in  two 
dimensions  by  NMDS  with  a  stress  value  =  0.01 
(Figure  SB).  Of  the  81  species  captured,  13  species 
were  captured  only  in  the  lower  pools  and  30  only 
in  the  upper  pools.  Nevertheless,  several  of  the 
species  captured  only  in  the  upper  pools,  including 
bigmouth  buffalo,  bowfin,  quillback,  and  walleye, 
are  known  to  be  established  in  the  lower  pools. 
Thus,  at  least  some  of  the  difference  in  community 
composition  between  upper  and  lower  pools 
reflects  sampling  artifacts. 
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Community  Structure  versus  Distance 


figure  6.  Relations  between  the  distance  between  pools  (river  miles)  with  (A)  community  structure  of  fish  (Bray-Curtis  similiarity),  (B) 
community  composition  of  fishes  (Euclidian  distance),  and  (C)  habitat  composition  (normalized  Euclidian  distance).  Plotted  are  all  possible 
pairwise  comparisons  among  the  six  regional  trend  areas  and  outpools. 


Figure  7.  Relation  between  water  depth,  flow,  conductivity,  vegetation  density,  presence  of  woody  structure,  and  presence  of  flooded  terrestrial 
vegetation  with  river  mile  (relative  to  Pool  3).  All  habitat  measures  were  collected  at  the  electrofishing  sites  sampled  in  2000. 
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UMRS  Navigation  Pool 


Figure  8.  W  Cluster  analysis 
of  fish  community  composition 
data  collected  through  seining. 
Major  pool  groupings  (upper  and 
lower)  are  noted.  Nonmetric 
multidimensional  scaling  plot 
and  stress  value  for  Upper 
Mississippi  River  System  pools  in 
two  dimensions.  The  two  major 
groups  from  the  cluster  analysis 
are  easily  identifiable. 


Discussion 

Our  analysis  of  both  the  community 
composition  and  community  structure  of  fishes 
in  the  UMRS  each  yielded  two  major  pool 
groups,  upjjer  and  lower.  Two  previous  studies 
also  classified  UMRS  pools  into  upper  and 
lower  reaches  based  on  habitat  variables  (U.S. 
Geological  Survey  1999;  Koel  2001).  It  is  likely 
that  geographic  range  limitations  of  fishes,  habitat 
factors,  and  possibly  historical  barriers  have  all 
influenced  the  fish  composition  and  community 
structure  differences  between  upper  and  lower 
pools.  Our  analysis  also  revealed  four  or  five 
subgroups  of  pools.  Based  on  the  strength  of  our 


analysis  of  similarity  tests  (i.e.,  Global  R  values) 
and  NMDS  plots,  these  subgroupings  may  present 
a  more  accurate  description  of  the  similarity  of 
community  composition  and  structure  among 
the  UMRS  pools  sampled.  Although  there  were 
clear  differences  between  the  upper  and  lower 
pool  groups  based  on  community  structure  and 
community  composition,  it  is  clear  that  more 
spatial  structure  exists  in  this  system  than  a  simple 
dichotomy  of  upper  and  lower  pools. 

Results  from  this  study  should  be  interpreted 
with  caution  because  the  data  available  for  analysis 
were  limited  to  1  year  (covariation  of  communities 
cannot  be  addressed)  and  essentially  one  sampling 
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gear,  an  electrofishing  boat.  Lubinski  et  al.  (2001) 
reported  that  among  all  gears  used  in  the  LTRMP, 
electrofishing  generally  had  the  greatest  statistical 
power  to  detect  trends  across  all  species  and  habitat 
types.  Nevertheless,  boat  electrofishing  does  not 
sample  all  species  within  the  UMRS  equally 
well.  For  example,  electrofishing  is  conducted 
near  the  shoreline  and  will  not  be  effective  for 
species  that  primarily  occupy  offshore  habitats. 
Because  of  the  limitations  of  boat  electrofishing, 
seining  was  included  in  the  plans  for  this  study  to 
provide  additional  information  on  the  small  fish 
community.  Unfortunately,  we  found  seining  to  be 
logistically  untenable  to  conduct  in  all  outpools, 
but  the  limited  seining  data  collected  showed 
major  pool  groupings  (upper  and  lower  pools) 
consistent  with  our  analysis  of  electrofishing  data. 
Given  the  vast  spatial  extent  of  the  UMRS  and  the 
great  diversity  of  habitat  types  and  fish  species  it 
contains,  it  is  unlikely  that  any  single  study  could 
fully  address  both  the  patterns  and  causes  of  spatial 
variation  of  fish  communities  within  this  system. 
We  feel  the  present  study  provides  a  useful  first 
step  in  addressing  this  issue. 

In  general,  outpools  tended  to  group  with 
adjacent  RTA  pools.  Outpools  5, 12, 14,  29, 
and  31  were  within  the  same  subgroups  as 
their  adjacent  RTA  pools  (4, 13,  and  the  Open 
River  Reach).  These  results,  and  our  habitat 
correlation  analysis,  suggest  a  strong  negative 
relation  between  the  distance  among  pools  and 
similarity  of  fish  community  composition  and 
structure.  In  other  words,  our  results  suggest  fish 
communities  in  adjacent  UMRS  pools  and  reaches 
tend  to  be  similar.  The  exception  to  this  trend, 
the  subgrouping  of  outpools  3, 19,  and  20  in  our 
community  stracture  analysis,  may  have  arisen  as  a 
result  of  similarity  in  habitat  features  or  low  sample 
sizes  in  Pools  19  and  20.  La  Grange  Pool  was  a 
unique  subgroup  for  both  community  composition 
and  structure  analyses,  which  was  an  expected 
result  for  this  tributary  RTA  pool.  Pool  8  was  a 
unique  subgroup  in  terms  of  community  structure. 
This  RTA  pool  had  the  greatest  abundance  of 
centrarchid  species,  which  may  be  related  to  the 
relatively  greater  abundance  of  aquatic  vegetation 
found  in  this  pool  (Figure  7). 


Our  attempts  to  correlate  spatial  variation  of 
fish  communities  with  habitat  data  were  hindered 
by  the  confounding  of  habitat  similarity  and 
distance  between  pools.  Both  the  composition 
and  community  structure  of  fishes  should  vary 
as  a  function  of  distance  between  pools  because 
of  zoogeography,  immigration  and  emigration, 
source-sink  dynamics  and  similar  histories  of 
large  scale  disturbances  such  as  major  floods  and 
droughts  (Drake  1990, 1991 ;  Hamrick  and  Nason 
1996;  Pullium  1996).  Because  habitat  similarity 
was  also  correlated  with  distance  between  pools, 
it  is  difficult  to  determine  the  influence  of  habitat 
on  fish  communities  independent  of  the  spatial 
demographic  processes  listed  above.  To  gain  a 
better  understanding  of  the  influence  of  habitat 
on  fish  communities,  future  studies  could  attempt 
to  account  for  both  spatial  proximity  and  habitat 
variation  by  selecting  pairs  of  study  pools  that  are 
relatively  close  together,  but  differ  substantially 
in  specific  habitat  measures.  Also,  future  analyses 
could  devise  an  index  of  historic  habitat  alterations 
for  each  RTA  and  outpool  to  assess  if  fish 
community  variation  correlates  with  this  index. 

This  study  was  not  able  to  resolve  where 
Navigation  Pools  19  and  20  fit  within  the  UMRS 
as  a  whole.  Pools  19  and  20  were  similar  to  lower 
pools  in  terms  of  community  composition,  but 
similar  to  upper  pools  with  regard  to  community 
structure.  Electrofishing  collections  from  Pools 
19  and  20  differed  from  all  other  pools  in  that  all 
data  were  collected  over  a  period  of  2  days,  rather 
than  over  a  period  of  5  months.  It  is  interesting, 
however,  that  Pools  19  and  20  grouped  together 
in  both  the  community  composition  and  structure 
analyses  because  Lock  and  Dam  19  is  known  to 
be  a  barrier  to  migratory  fishes  such  as  skipjack 
herring  (Kelner  and  Seitman  2000).  Despite 
this  barrier,  this  study  suggests  the  overall  fish 
communities  with  Pools  19  and  20  are  relatively 
similar.  An  important  caveat  to  this  study  is  that 
only  three  UMRS  pools  below  Pool  14  were 
sampled  (19, 20,  and  26),  whereas  seven  upper 
UMRS  pools  were  sampled  (3, 4,  5,  8, 12, 13,  and 
14).  Studies  including  a  greater  number  of  lower 
UMRS  pools  might  improve  our  understanding  of 
spatial  variation  of  fish  communities. 
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Implications  for  Long  Term  Resource 
Monitoring  Program  and  Future  Studies 

Current  RTA  areas  are  evenly  distributed  within 
the  major  pool  groupings  identified  in  this  study 
(i.e.,  three  RTA  areas  in  upper  pool  group,  three 
RTA  areas  in  the  lower  group),  which  supports  the 
premise  that  LTRMP  fisheries  data  can  be  used  to 
make  inferences  to  the  entire  UMRS.  Subgroupings 
of  outpools  with  nearby  RTA  pools  and  the 
importance  of  distance  between  pools  in  habitat 
correlations  suggest  that  fish  community  data  from 
RTA  pools  should  at  least  be  relevant  to  other 
nearby  UMRS  pools.  Furthermore,  these  results 
suggest  that  expanding  LTRMP  fish  monitoring  to 
pools  adjacent  to  current  RTA  areas  would  yield 
minimal  additional  information.  Further  research 
is  needed  to  resolve  how  fish  communities  within 
Navigation  Pools  19  and  20  and  other  lower  UMRS 
pools  compare  to  current  RTA  areas.  Future  studies 
in  this  area  should  (1)  address  covariation  of 
community  measures  through  time,  (2)  examine 
additional  pools  and  reaches  in  the  lower  portion 
of  the  system,  and  (3)  further  examine  relations 
between  fish  communities  and  habitat  features 
using  experimental  designs  that  specifically 
account  for  the  confounding  of  habitat  similarity 
and  distance  between  pools. 
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